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j3 ' Abstract. We discuss how the spectral changes of quarkonia at T c can reflect 

the "critical" behaviour of QCD phase transition. Starting from the temperature 
dependencies of the energy density and pressure from lattice QCD calculation, we 
extract the temperature dependencies of the scalar and spin 2 gluon condensates near 
T c . We also parameterize these changes into the electric and magnetic condensate near 
T c . While the magnetic condensate hardly changes across T c , we find that the electric 
condensate increases abruptly above T c . Similar abrupt change is also seen in the 
scalar condensate. Using the QCD second-order Stark effect and QCD sum rules, we 
^ (— 1 show that these sudden changes induce equally abrupt changes in the mass and width 

of J/ip, both of which are larger than 100 MeV at slightly above T c . 



1. Introduction 

O. 

I/-} . J /if> emanating from heavy ion collision is considered to be an intriguing tool to 
investigate the properties of quark gluon plasma (QGP) expected to form in the early 
stages of a heavy ion collision (TJ [2] . Indeed measurements at past SPS and recent RHIC 
data show nontrivial suppression patterns that could be consistent with the original 
^ . prediction that the J/ip will dissolve in QGP. However, the subject has gained a new turn 
recently as lattice calculations show that J / ip will survive past the critical temperature 
T c up to about 1.6T C , while \ c and ip' will dissolve just above T c 0, Hj. The notion 
that Jjijj will not dissolve immediately above T c was also considered before [5], as it 
was known that non-perturbative aspects of QCD persist in QGP [SI EJ- These results 
suggest that the sudden disappearance of J / ip at T c is not the order parameter of QCD 
phase transition. On the other hand, the phase transition is accompanied by sudden 
changes in the chiral condensate, the heavy quark potential and the energy density, and 
should therefore have some effect on the properties of quarkonia at T c . 

Unfortunately, the lattice results based on the maximum entropy method have large 
error bars, and one can not determine the detailed properties of quarkonium above T c , 
except its existence. Solving the Schrodinger equation with potential extracted from the 
lattice is one possibility |8]. However, here there are some controversies on how to extract 
the correct potential from the free energy calculations. Even with a chosen prescription, 
one has to solve the Schrodinger equation at discrete temperatures steps, and it becomes 
not clear how the sudden critical behaviour is translated into the discrete temperature 
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steps. Thermal perturbation alone is also not sufficient to probe the temperature region 
from T c to 2.5T c [9j, which is now known to be strongly interacting. 

Therefore, a systematic non-perturbative method is essential to treat the phase 
transition region. Here, we will summarize the recent developments [TOl [TT] to attach 
the problem using QCD sum rules and the QCD second order Stark effect. The inputs 
are the temperature dependencies of local gluonic operators, which undergoes abrupt 
changes across the phase transition as does the energy density. 



2. The gluon matter 

The order parameter for QCD phase transition is the thermal Wilson line for pure 
gauge theory, and the quark condensates for QCD with massless quarks. For realistic 
QCD none of them are order parameters in the strict sense. But simulation shows that 
the sudden changes in both parameters take places at the same critical temperature 
T c , which is determined by the susceptibilities. At the same T c , the energy density 
also makes a drastic change, which is somehow universal for any quark flavours. The 
sudden changes apparent in the energy density and pressure can be translated to 
temperature dependencies of local gluonic operators, which are expected to embody the 
non perturbative nature of QCD. The link is obtained through the energy momentum 
tensor, which can be written in terms of symmetric traceless part (twist-2 gluon) and 
the trace part (gluon condensate). 



ST(G aaX Gf ) + ^^-G%G a ^ } (1) 



where we will use the LO beta function (3(g) = — — §iVf). The thermal 

expectation value of this operator can be related to the lattice measurements of energy- 
momentum tensor through the following relation. 

(T^) T = (e + p) fuV - \gA + - A (e - 3p)g a ^. (2) 

Here, u a is the four velocity of the heat bath. With the following definitions, 

S -G%G a ^\ = G (T) (3) 

! ST(^G aaX Gf)) T = (V?/ - G 2 (T), (4) 

we find for pure SU(3) gauge theory, 

G (T) = G™-^-(e-3p), G 2 (T) = -^-(e + p), (5) 

11 7T 

where G™ c is the value of the scalar gluon condensate in vacuum. 

At low temperature, the scalar part is dominated and characterized by the non- 
perturbative contribution of the QCD vacuum. At high temperature, its behaviour will 
be dominated by the quasi-particle contribution. Just above T c , it is still dominated 
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Figure 1. a) Gluon condensates and b) electric and magnetic condensates near T c . 



by the non-perturbative contribution. For the heat bath at rest, one can rewrite the 
thermal expectation values in terms of electric and magnetic condensate. 

(t^t = - 3 G " (T) " I G!(T) (6) 

(V fi2 )r = \ G " iT) - i G2(T) ' (7) 
Figure 1 shows the temperature dependence of Go, G-i or E 2 and B 2 . One should first 
note that the sudden increase of energy density is translated to the anomalously large 
and sudden decrease in Go, which deviates largely from the asymptotic T 4 behaviour. 
In contrast, G2 reaches the asymptotic value quickly Similarly, one finds that there 
is a sudden increase in the electric condensate E 2 , while the magnetic condensate B 2 
hardly changes above T c . This can be related to the fact that the area law behaviour 
of the space time Wilson loop changes to the perimeter law above T c , while that of the 
space-space Wilson loop does not [6]. The connection comes in through the operator 
product expansion (OPE) of rectangular Wilson loop, which was found to be expressible 
in terms of the electric condensates and the magnetic condensates for the space-time 
and space-space Wilson loops respectively [12]. To investigate the consequences of the 
abrupt changes of condensates to the properties of J/i[>, we will use perturbative QCD 
and QCD sum rules. 

3. Results from perturbative QCD 

The perturbative QCD formalism for calculating the interaction between heavy 
quarkonium and partons was first developed by Peskin [131 EJ- The counting scheme 
in this formalism is obtained in the non-relativistic limit. Because of this, the formula 
for the mass shift reduces to the second-order Stark effect in QCD, which was used to 
calculate the mass shift of charmonium in nuclear matter [15]. The information needed 
from the medium is the electric field square. Noting that the dominant change across the 
phase transition is that of the electric condensate, one finds that the QCD second-order 
Stark effect is the most natural formula to be used across the phase transition. 
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Figure 2. Mass shift from the second-order Stark effect (solid line) and the maximal 
mass shift obtained from QCD sum rules from Ref. [10] (points). 



The QCD second-order Stark effect for the ground state charmonium with 
momentum space wave function normalized as J ?f^3" IV^P)! 2 = 1 is as follows, 
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(8) 



where k = \k\ and (^AE 2 )? denotes the value of change of the electric condensate 
from its vacuum value so that 5m j« = in vacuum. The second line is obtained 
for the Coulomb wave function. Here, e is the binding energy and m c the charm quark 
mass. These parameters are fit to the size of the wave function obtained in the potential 
model [16], and to the mass of J/ip assuming it to be a Coulombic bound state in the 
heavy quark limit [13J. The fit gives m c = 1704 MeV, a = 0.271 fm and a s = 0.57. 
Few comments are in order. The minus sign in (jSJ) is a model independent result and 
follows from the fact that the second-order Stark effect is negative for the ground state. 
Second, the important parts of the formula are the Bohr radius a, the binding energy e 
and the temperature dependence of the electric condensate squared. In the formula, the 
factor of a 2 follows from the dipole nature of the interaction, and the binding energy 
from the inverse propagator, characterizing the separation scale [T3], [T7]. Therefore, the 
actual value of the mass shift does not depend much on the form of the wave function 
as long as the size of the wave function is fixed. 

The solid line in Fig. [2] shows the mass shift obtained from the second-order Stark 
effect with the extracted lattice input for the electric condensate shown in Fig. 1-b. As 
the formula is based on multipole expansion, it will break down if the higher dimensional 
condensates become large. At this moment, there is no lattice calculation for the 
temperature dependence of higher dimensional gluonic operators. However, as discussed 
before, up to 1.1T C the temperature dependence is dominated by the sudden increase in 
the electric condensate with little change in the magnetic counterpart. Moreover, the 
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OPE in the QCD sum rules were found to be reliable up to 1.05T C [10]. Therefore our 
result should at least be valid up to the same temperature. Above that temperature, 
the change of the electric condensate amounts to more than 100% of its vacuum value 
and higher order corrections should be important. As can be seen in Fig. 1, the mass 
reduces abruptly above T c and becomes smaller by about 100 MeV at 1.05T C , reflecting 
the critical behaviour of the QCD phase transition. 

4. QCD sum rule result 

The QCD sum rule for quarkonium was found to be very reliable at zero temperature. 
This is because the expansion parameter of the OPE for the correlation function for the 
heavy quark current-current correlation function is Aq CD /(Q 2 + Am 2 ). Therefore, the 
OPE for the correlation function can be reliably estimated even at Q 2 = 0. Moreover, 
for the heavy quark system, only gluon operators appear, for which reliable estimates 
can be made at least for the lowest dimension 4 operator. 

How to generalize the sum rule method to finite temperature depends on the 
magnitude of the temperature. If the temperature is larger than the separation scale 
Q 2 + 4m 2 , we have to calculate the temperature effect into the Wilson coefficient. If the 
temperature is low, then all the temperature effects can be put into the temperature 
dependent operators. In this case, the new expansion parameter in the OPE will be 
(Aqcd + cT) 2 /(Q 2 + Am 2 ), where c is some constant. Whether such approximation 
is valid or not can be directly checked by looking at the convergence of the OPE at 
finite temperature. An additional ingredient at finite temperature is that unlike at zero 
temperature, where only scalar gluon operators appear, operators with Lorentz indexes 
have to be added. Therefore, up to dimension 4, both the scalar gluon (Go) and twist-2 
gluon {Gq) operators contribute. Since we have extracted the temperature dependence 
of these operators from the lattice, no ambiguities exist for the OPE side of the sum 
rule for a consistent analysis. 

The OPE for the correlation function of heavy quark currents is related to the 
phenomenological side via the dispersion relation. In the sum rule, one has to assume 
the form of the spectral density. Here, one notes that lattice calculations suggest that the 
peak structure in the spectral density persists above T c , although the resolution is not 
good enough for a critical study. Therefore, we can assume a relativistic Breit-Wigner 
form for the spectral density. Using the moments for the correlation function, we can 
extract a constraint for the temperature dependence for the mass and width from the 
temperature dependence of the condensates. We find that if there is no change in the 
width, the mass critically decreases by few hundreds MeV slightly above T c , equivalently 
if there is no change in mass, the width will increase by a similar amount pU]. The dots 
in Fig. 2 represent the maximum mass shift assuming no change in the width. This is a 
direct consequence of the critical change of the scalar condensate proportional to e — 3p. 
As can be seen in the figure, the mass shift obtained from the second-order Stark effect 
is almost the same as the maximum mass shift obtained in the sum rule up to T c and 



Properties of Quarkonia at T c 



6 



250 

200 

? 150 

a> 

S 

U 100 
50 




<Xqq> 5=1 • 
a qq > 5=2 
• a aa >5=3 ♦ 




4 






♦ 














JA|/ 
♦ ■ 
• 


i 






• 






ft 







0.9 



0.95 



1.05 



1.1 



T/T„ 



Figure 3. Thermal width of J/-0 obtained from the second-order Stark effect and 
QCD sum rule constraint. 



then becomes smaller. The mass shift at T c is about —50 MeV. In the QCD sum rules, 
only a constraint for the combined mass shift and thermal width could be obtained. 
This constraint can be crudely summarized as follows, 

- Am + T T ~ 80 + 17 x (T-T c ) [MeV], (9) 

within the temperature range from T c to 1.05T C . Therefore, the difference between the 
Stark effect and the maximum mass shift obtained from QCD sum rules above T c in 
Fig. [2] could be attributed to the non-perturbative thermal width at finite temperature. 
In Fig. El we plot the thermal width obtained from combining the constraint in (jHJ) 
with the mass shift obtained from the QCD second-order Stark effect. As can be seen 
in the figure, the thermal width at 1.05T C becomes larger than 100 MeV. Such width 
slightly above T c is larger than that estimated from a perturbative LO and NLO QCD 
matrix element calculation together with an assumption of weakly interacting quarks 
and gluons with thermal masses jT8j [19], but smaller than a recent phenomeno logical 
estimate [20] . 

The mass of quarkonium at finite temperature was also investigated in the potential 
models [2T] . where the mass was found to decrease at high temperature. However, 
the potential has to be extracted from the lattice at each temperature and hence a 
more detailed investigation are needed to identify the critical behaviour of J/ip in the 
temperature region near T c . 

5. Corrections from dynamical quarks 

To consider the quark effects, first we consider the quark operators appearing in the 
OPE. We can neglect the light quark contribution to the OPE, because the light 
quark operators appear in the OPE at order a 2 s (q 2 ): This is why the light quark 
condensate can be neglected in the sum rules for heavy quark system. On the other 
hand, thermal heavy quarks that directly couple to the heavy quark current contribute 
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to the OPE at leading order. This is different from the heavy quark condensates that 
are perturbatively generated from the gluon condensates, and contribute to the OPE 
through gluon condensates, whose Wilson coefficients are calculated in the momentum 
representation. The direct thermal quark contributions are called the scattering terms. 
However, similar terms also appear in the phenomenological spectral density, which also 
has free charm quark mode that is not coupled with a light quark in the form of a D 
meson above T c as been recently studied in Ref. [22]. Therefore, the scattering term will 
cancel out between the OPE and the phenomenological spectral density in the QCD 
sum rule analysis for the deconfined medium. 

Second, the gluon condensates itself can have a different temperature dependence 
in the presence of dynamical quarks. As discussed before, the important input for the 
mass and width change is the temperature dependence of gluon condensates in Fig. [TJ 
in particular the dominant contribution comes from the temperature dependence of Go- 
For that purpose, we note that the trace of the energy momentum tensor to leading 
order is given as, 



Therefore, we start from the lattice calculation of the trace of the energy momentum 
tensor for the full QCD with realistic quark masses as given in Ref. |23j. Then, we 
subtract the fermionic part of the trace anomaly, which was also shown in the literature, 
from the total. Then we divide the result for the relevant prefactor with Nf — 3 
multiplying the gluon condensate as given in ( flOl) . Since the critical temperature T c 
differs, we compared it as a function of T/T c in which T c = 196 MeV for the full QCD 
case [23J. As can be seen in Fig.HJ the magnitude of the resulting change near the critical 
temperature are remarkably similar between the full and pure gluon QCD; although the 
slope at T c is milder for full QCD as a consequence of rapid cross over transition instead 
of a first order phase transition. Since the change of the condensate sets in at a lower 
T/T c in the full QCD case, the mass and width of charmonia might start varying at 
a lower temperature in the realistic case than in the pure glue theory. Therefore we 
believe our main argument and the quantitative result will not be alter even in the 
realistic situation. 

6. Summary 

We have shown that the critical mass shift and width increase of quarkonia slightly above 
T c could effectively be an "order parameter" of QCD phase transition. The expected 
shift for J /if) is small and will be a challenge for LHC. However, larger shift is expected 
for Xc at its formation temperature slightly above T c . Therefore, direct measurements 
and confirmation is possible. The changes will also lead to changes in production ratios 
in the statistical model and J /if) suppression effects, which needs further detailed studies. 
This work was supported by BK21 Program of the Korean Ministry of Education, and 
by the Korean Research Foundation KRF-2006-C00011. 
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Figure 4. Comparison of the scalar gluon condensate in the pure gauge theory with 
the one extracted from full lattice QCD. Horizontal errorbars in the full QCD case 
are drawn by assuming the 2% uncertainty in the conversion from the lattice units to 
physical temperature [23) . 
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